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We introduce the chemical potential in a system of two-flavored massless 
fermions in a chiral bag by imposing boundary conditions in the Euclidean 
time direction. We express the fermionic mean number in terms of a functional 
trace involving the Green function of the boundary value problem, which 
is studied analytically. Numerical evaluations for the fermionic number are 
presented. 



OO 
ON 




12.39.Ba, 11.10. Wx, 12.38.Mh 




1 



Typeset using REVTpjX 



I. INTRODUCTION 



Functional determinants of elliptic differential operators are useful for describing one-loop 
effects in Quantum Field Theory and Statistical Mechanics. 

When the fields are constrained to live inside a bounded manifold in the Euclidean 
space-time, boundary conditions become relevant in the calculation of physical quantities. 

This is the case of Hybrid Chiral bags, which are effective models for confinement [l|-|3| 
(For a general review on bag models see M). They share the characteristics of two succesful 
approaches to describe baryons: bag models |||6| and Skyrme model [0-f|. 

According to the Cheshire Cat Principle (CCP) PIJTH, in these two-phase models, 
fermionic degrees of freedom could be replaced by bosonic ones in a region of space, physical 
quantities being insensitive to the position of the limit of separation of the two phases. In 
I + I-dimensions the CCP follows from the bosonization of fermionic fields JT]|. In the 3 + 1- 
case, topological quantities such as the baryonic number at T = have a similar behavior 



11], while for non-topological ones, as the total energy at T = 0, the CCP is approximately 



valid @. 

In this paper we will study the mean fermionic number for a two phase model at finite 
temperature and chemical potential. This represents an extension of previous calculations 
done in the frame of the MIT bag model . 

In Section we will set the problem, expressing the mean fermionic bag number as an 
adequate trace of the Green function of the system. In doing this, the chemical potential is 
introduced through the boundary conditions in the Euclidean time direction |14|,[13|. 

Section [III] is devoted to the calculation of the Green function for the Dirac operator 
obeying the above temporal boundary condition, and spatial (chiral bag) boundary condi- 
tions which relate the fermionic degrees of freedom in the inner phase with the Skyrme field 
in a hedgehog configuration outside of the bag. 

In Section |V| the expression of the mean fermionic number is built. By means of the 
asymptotic Debye expansion of Bessel functions we are able to split the mean fermionic 
number into two pieces: a regular one, to be evaluated numerically, and a possibly non- 
regular one, to be treated analytically following the techniques presented in Ref. [|T3] . 

In Section [V| we show our numerical results, for the total mean fermionic number includ- 
ing both the contributions of the bag sector and of the skyrmion sector as in Ref. M. 

Finally, in Section [Vl the conclusions are presented. 



II. THE HYBRID BAG MODEL 

We will be interested in the study of the mean fermionic number for a two phase model: 
a two-flavored massless fermionic field confined inside a static sphere of radius R and a 
hedgehog Skyrme field filling the exterior sector. These two phases are linked by the spatial 
chiral boundary conditions, which will be introduced in the following. 

The fermionic field inside the bag S <g> [0,1], < \x\,t < 1 is described by the Dirac 
operator 



D(fi,R) 



It 



(1) 
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Here X/ denotes the identity in the flavor (isospin) space, and (3 stands for the inverse 
temperature. 

The external phase is described by a Skyrme model M,T§\ , whose Lagrangian is given by 



C 



: F^Tr (d^Ud^U^ + -^Tr 



16 n 



32e 2 



(d.U) U\{d v U) 



(2) 



where the scalar field U(x) takes values in the SU (2) group. 

When the whole space is filled by U(x), a topological stable classical solution is given by 
the hedgehog configuration 



Un 



(3) 



where the profile of the skyrmion 9(r) is the chiral angle. In this pure skyrmionic model, 
the imposition of the boundary conditions 



0(r = 0) = 7T, 9{r 



oo 







leads to a topological charge (winding number) B — 1, which is identified with the baryonic 
number of the skyrmion The profile of the chiral angle, 9(r), is numerically determined 
by minimizing the classical energy of the soliton P,|l6 . 



When the skyrmion fills only the exterior of the bag, its contribution to the fermionic 
number, iVsk, is given by the integral of the topological charge density over that region, 
giving 

7T 



Sk 



sin 9 cos ( 



(4) 



as discussed in Ref. [||. 

In the hybrid bag model, fermionic and skyrmionic phases are connected through the 
spatial boundary conditions 



Bip(t,x) = for \x\ = 1 



where 



B = X 4 ® X/ + (z jk®Ti) e 



-i97 5 ®(r-n) 



(5) 



(6) 



being 9 = 9(R) the value of the skyrmion profile at the surface of the bag. 

Since we are dealing with fermions, these fields must be antiperiodic in the Euclidean 
time direction, 



^(1, x) = — ^(0, x). 
The Grand Canonical partition function is given by 

S(T, R, fj.) = e -P G{T ' R ^ ~ Det \D((3, R) - t^ 



BC 



(7) 



(8) 



where "BC" means that the differential operator is defined on a space of functions satisfying 
both spatial and temporal boundary conditions. In the above expressions, fi stands for the 
chemical potential. 
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As discussed in Ref. [13| the mean fermionic number can be expressed as 



BC 1 

(N hag )((3, R, li; 9) = -—(/?, R, li; 9) = -Tr{-i ( 7 ° ® X/) k(t, x; t', x')}, 

where k(t, x; t', x') is the Green function of the problem 

D(f3,R)k(t,x;t',x') = 8^(x-x')S(t-t')I 8 , forx,x' G E, 
Bk(t, x; t', x') = 0, forx G <9E, 

fc(l, f ; t', x>) + e _/l/? fc(0, f ; t' , x>) = 0. 

The total mean fermionic number for the hybrid chiral bag model is defined as 

(N) (/?, R, li; 9) = (iV bag ) (/?, R, li; 0) + N Sk (9) 



(9) 



(10) 

(11) 
(12) 



(13) 



One can apply discrete transformations to these equations to obtain symmetries of (N^ ag ). 
In fact, 

ki(t, x; t', a?) = -7 5 A;(t, x; t', x')rf (14) 

is the Green function of the problem with 9 — > 9 + tc in the spatial boundary condition, 
equation ([TT|) . Then, from equation (||), it is easy to establish that 

(iV bag ) (ft fl, + 7T) = (iV bag ) (/?, i2, ^ 0) (15) 
Taking into account (|4]) and (IS), the total mean fermionic number satisfies 

(N) (p, R,LL-9 + n) = l + (N) (p, R, la; 9). (16) 

Moreover, 

k 2 (t, x- 1', x>) = 7° 7 5 fc(l - t, x ; 1 - x'hV (17) 

is the Green function of the problem with 9 — ► — in the spatial boundary condition and 
/Lt — > — /i in the Euclidean time boundary condition (|T2"D . So, from @ it can be seen that 

(N hag )(P,R,-fi;-9) = -(N hag )(f3,R,Li;9). (18) 

Notice that the total mean fermionic number has the same behavior, since Nsk(—9) = 
-N Sk (9). 

From the above relations, one can deduce that 

(N)(P, R, -li; 9) = l- (N)(P, R,Kn-0), (19) 
allowing to determine (JV) for li < from its values for \i > and < 9 < n. 
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III. GREEN FUNCTION 



In this section we will obtain the Green function in equations (|I0),(O),(^). In order to 
satisfy the "temporal" boundary conditions, we propose for the Green function 

oo 

k(t,x;t',x') = R M£ > ^)e~ <nn/3(t " t ' ) , (20) 

n=— oo 

where 

7T 

fl n = u n -ip, u n = (2n + l)-. 

We adopt the representation for the Euclidean Dirac matrices, 

7° = ip 3 ®X 2 ^ = ip 2 ®a 7 5 = p 1 ®J 2 , (21) 
where a h , p k are Pauli matrices. For each value of n, we obtain 

iQ n p 3 <g> J 2 <g> Zi - p 2 ® a ■ V ® X/] fc n (f , x') = <5 (3) (f - a?) J 2 ® X 2 ® 2}, (22) 

with f2 n = f2 n i? = (2n + 1)712; — z/i, being 2; = i?T and p = pR dimensionless variables to be 
used along the paper. 

We write the Green function as the sum of a particular solution of the inhomogeneous 
differential equation plus a solution of the homogeneous equation, 

') = k^(x,x') + ~k (x,x'), (23) 

with 

3 

fc£» (x , x>) =Y.P l ® Ai£,n) & , (24) 
3 

k n (x, x>) = <8> a {e ' n) (x, x>). (25) 
e=o 

We choose k$(x, s') so as to give the contribution of a free gas to the mean fermionic 
number. The second term (p5|) accounts for the correction, due to the boundary conditions, 
to the mean fermionic number. 

The differential operator ([!]) and the boundary condition operator (g|) are invariant under 
transformations in the diagonal subgroup of SU{2) rotations ® SU(2) isoS p in . So, the Green 
function is block-diagonal and can be expressed in each invariant subspace in terms of the 
eigenstates of {K 2 , J 2 , L 2 , S 2 , 1 2 , K 3 }, where K = L + S + I. For this purpose, we choose 
the following basis in the (k, m) subspace: 

|1) = \k, J = k + 2, 1 = k, m) 1 2) = \k, j = k — |, Z = k,m) 
|3) = \k,j = k + l,l = k + l,m) |4) = \k,j = k - f, I = k - l,m) 
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for k 7^ 0, and 

\l)o = \0,j = U = 0,0) 
|3) = |0,j = i,/ = l,0) 

for k = 0. 

The free field contribution to the Green function in the invariant subspaces labeled by 
n, k, m is given by 

A$ m) (ry)=A$ m) (r,r') = (26) 
A< (k,m)(^ r ') =iSn&lJ(k, m )(r < )H {k , m) {r > ) (27) 



( i \ _q2 J ■/(fc,m)(r-<)M( fc)m )^ (fc>m )(r>) r' > r , , 

( fc ' m ) 1 ' J " n \ ff( fc , TO )(r>)M (fc>TO) J (fc>m) (r<) r' < r ' ^ 

where the dependence on n is implicit in the arguments of the spherical Bessel functions 
appearing in the expressions for the diagonal matrices 

J(k,m)(r) = diag (j k (iS n (i n r), j fc (^Ar),jVi(*£A r )» jk-i(iS n Q n r)^ (29) 

J(Q,0)( r ) = dia § (MiSn&nrfJliiSnClnrf) (30) 

H(k, m ){r) = diag (hljp (iS n Q n r), h[ 1] \iS n Q n r), h^l^iSAr), h^l^iSAr)) (31) 

H m (r) = diag (hS\iS n ti n r), h?(tS n Q n r)) . (32) 

In previous equations, S n stands for the sign of u> n and 

M {Km) = I ° ^ ) , M (0 ,o) = ( I I J . 



We propose the following Ansatz for the coefficients introduced in fl25|) when considering 
the solution of the homogeneous differential equation, 

a (°!m)( r > r/ ) = iSn&lJ {k , m )(r)c(0)^ m) J {k , m) (r') 

a< (k]m)( r ^ r ') = -^n^(fc,m)( r ) M (fc, m )C(0)( fcim )J( fc>m )(r / ) 
a (fc!m)( r ' r ') = -^n^(fc,m)( r ) M (fc,m)C(3)( fcim )J( fcim) (r / ) 

a (fc!m)( r ' r ') = iS n £llJ( kim )(r)c(3)( k , m ) J( fc , m )(r'). (33) 

By imposing the boundary conditions, we can determine the unknown matrices c(0V& m ) 
and c(3)(fc )7n ) from the equations 

i( 2 ) ™a i .( 2 ) rn i i^aast z^t\ „(°) 



A \k,m) & x ') + G (fc,m) & x ') + icos6{in-a® Ti) (k>m) a\ k \ m) (x, 

) a {1) 



-ism 6 {in-a®in-r) (k m] a ( y m] {x,x') = 0, (34) 
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A { n}^ (x, x') + ag m) (x , x') + cos 9 (in-a®li) (k ^ (x , a?) 



- sin {in- a® in- f) (fem) a[°? m) (f, 2?) = 0, (35) 



for |x| = 1. The matrices (in • a (g> Ij) * and (in ■ a ® in ■ T) (km) are explicitly given in 



Ref. g7|. 

It is straightforward to see that, for k ^ 0, 



-(0)(fc,m) 



^0 U 2x2 
^2x2 ~~ 



U9x9 — L/ n 



"(3)(fc,m) 



(fe,m) 



cj; m) o 

02x2 C*3, 



2x2 



(fc,m) I > 



where 



„ . _ sin # a . . _ . 

5 n cos 6/ J 2 H + — srn 9 E k A k 

v 2v 



M k ' m ) — n 



(36) 



(37) 



°3,± 



_i sin# a _, 
5 n cos 6 1 0-3 J ± J T ± — — a 3 + — sin 9 J, ai J± 
2i/ 2z/ 



o ? 



(38) 



Here we have introduced a = ^Au 2 — 1, z/ = + 1/2, and 



e fc \ 
-e fc _! j ' 



4 \ 
j' 



Q J _ jfc-l(»Sn^n) 

ifc + l(«5'n^ II ) 

Jfc-l(*Snfln) / 



with 



j fc+ i(z5 , w Q w )j fc (z5' w n n ) 
Jfc+i (iS n £l n ) — jf,(iS n VL n ) 

Moreover, for = 0, 



jk+i(iS n fl n )h k \iS n Q n ) — jk^SnQ^h^^iSnfln) 



(1) 



4 



/ (7(0.0) Q \ 
c(0) (0 ,o) = I ° Q _ C r(o,o) J > c ( 3 )(o,o) = 

and the coefficients are determined by the equations 

(S n cos9 +2sin^e )C ( S ' 0) = d, 



C<°; 0) 



o a 



(0,0) > 



(39) 



^(0,0) 
°3,+ 



S n cos0 (io^n^n))" 1 ^^^) + sin ^°' 0) - (j Q (iS n Cl n ))- 1 h { 1 > (iS n Cl n ), (40) 



(0,0) 



S n cos9(j l (iS n ti n ))- 1 jo(iS n ti n )-sm9 - ^(z^))- 1 ^^^)- (41) 



-1 /.(I)/ 



The above equations allow for the complete determination of the Green function. In 
the next section, we will give the explicit expressions for those coefficients needed for the 
evaluation of the mean fermionic number of the bag. 
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IV. THE MEAN FERMIONIC NUMBER 



The mean fermionic number for the bag in the hybrid model can be written as 

(N hag ) = {N ) + {N). (42) 

In equation fl42]) , (N ) is the contribution corresponding to a free fermionic field. It has 
been calculated in fllB"], to give 



47T 



(N Q ) = N f x- Lu^ + Ei 
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7T Z 



(43) 



with the number of flavors Nf = 2 in our case. See Ref. for an evaluation of (iVo) 
employing techniques similar to the present paper. 



Finally, we will evaluate the second term of (^) as 



W = ft 



-i (7 ® J/) k(t,x; t\x') 



(44) 



It is seen from the Dirac matrices representation (pl|) that only a^(t, x; t', x') must be 
considered. 

In taking the trace, the integration over angular variables can be done to obtain 



(N) = 2z r 2 dr dtt Y tr{a^ n) (x,x')} 
Jo Jn , n __ 



00 „i 
2z V iS n tt 2 n / r 2 dr 



n=— 00 
00 A; 



tr (■/(o,o)(r)c(3)(o j o)J(o,o)('" / 



+ E E tr (j(k, m )(r)c(3)^ m )J {Km) (r'] 



k=l m=—k 



(45) 



It is easily seen from Equation (45|) and from the definitions (^,^) that only the diagonal 



a 



elements of c(3)( k , m ), 
equations ([37l -[4l"l) leads to the equalities 



°3,± 



(k,m) 
3,± 



C^°f^ are needed. Their explicit evaluation from 



V,(A:,m)~ 
°3,+ 


11 


'Q(k,m)~ 


11 ' 


V,(fc,m)" 
°3,+ 


22 


'£t{k,m)~ 



22 



^(0,0) _ ^(fc,m) 
°3,+ — °3,- 5 



(46) 



allowing to write for (N) 

CO 

(N) = 2z J2 iS n& n 

n=— 00 
00 fc 

+ E E {[C: 

k=l m=—k 

+ [a 



dr [C^f j Q (iS n tt n r)j (iS n tt n r') + ji(z5' n Q n r)j 1 (z5' n ^ n r / ) 



(k,m) 
3,+ 



11 



j k (iS n tt n r)j k (iS n tt n r') + j k+ i{iS n Q n r)j k+1 (iS n tt n r') 



i(k,m) 
'3,+ 



j k (iS n tt n r)j k (iS n tt n r') + (^-iW) , • (47) 



S 



Notice that we have changed the order in which the series and the integral are taken. To 
properly define the resulting series we will take r' = r(l — e), keeping e > up to the end of 
the calculation. 



In equation (f|7]), we can do the m-sum, to get a factor 2v. We can also rearrange the 
/c-sum of the first term to include the k = subspace contribution. In this way, we obtain 



oo oo 



n=-oo k=Q 

oo oo 

lim £ 7rzCl n £ 

e ~ * n=— oo fc=0 



a 



11 



+ 2(u + l) 



a 



(fc+l,m) 



22 



(48) 



where we have introduced 

I„ = (1 - e)~ 1/2 C rdrJ v (iSn£l n r)J u (iS n £l n r(l - e)), 



with J^(x) the cilindrical Bessel functions. 
One can see that 



I„ + I 



(1 - e)- 3 / 2 J„(iSA(l - e))J u (iSA) r 



1/ T J-iz+1 



^ - dp + e Or 



— V 



(49) 



where = tS n Q n (l - e) [£ log J u (x)^ , s () ; ) 

The first factor in the fc-series on the r.h.s. of (f48|) contains all the dependence of (N) 
on the chiral angle 8(R). Using the recurrence relations for the Bessel functions, we can 
express this factor in terms of modified Bessel functions, 



2v 



°3,+ 



11 



+ 2(z/ + 1) 



a 



(fc+l,m) 
3,+ 



22 



2i\ x 2 
it) (J v (iS n Ct n )y 



_ M„N V + P V P V 
2v — -. — — h 



-2(u + l 



M v+l N v+l + P 2 +1 



where 



= «Sra COS [it — U + w] 



N u = S n cos 9 \—u — v — w] 



ism6 VT^¥ 



2v 



v - 2t 2 w 



i sm 



2u 



v + 2t 2 w 



M y = S n cos 9 



i sin 9 Vl^t 2 
t 



V4z/ 2 - 1, 



= cos 9 



- Or + 



t 2 
i-t 2 



v) — V 



t 2 



-isin9 Vl t2 (di=° + 



P v = -i sin 9 ^ 1 ^ dT°y/Au 2 - 1. 



(50) 

(51) 
(52) 
(53) 
(54) 

(55) 

(56) 
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Here, 



v 



p — v ~\~ x , t — , x — S n Q n , 

P 



u = r v (s n Q n )K' u (s n n n ), v = - (i„(s n n n )K v (s n n n )) , w = ^i u (s n n n )K v (s n n n ) 

x K ' x z 



The double series ([48|) is not absolutely convergent for e = 0, so we must keep e > 
in the general term up to the end of the calculations. A way to isolate this non regular 



behavior consists in the subtraction of the M-order asymptotic (Debye) expansion JT9| of 
the general term Nk >m (z, fi, 6; e), Af f m , as discussed in Ref. |[3 . 
So, we can write 



(N) =N 1 + N 2 

oo oo 



E E [Nk, m {z, Ml £ = o) - A{ m (z, M;e = 0) 



(57) 



n=— oo k=0 

oo oo 



+Sg E E A £m(*.M;e) 



n=— oo fc=0 



For M > 3 the e — > limit can be taken inside the double sum for N\ . In equation (p7\ ) 
we have used the expansion up to the order M = 6 to improve the numerical calculations to 
be done later. The non regular behavior of the series for (N) is isolated in the second term 
N 2 , which can be studied analytically. 

From equations (^9| , [50|) , by defining 



Fi(z,fi; e) = D.D.® 



Iv(S n ti n (l -e)) 
I u (S n £l n ) 



(58) 



and 

F 2 (z,p,,9;e) 



(2iz)Q n 



'1 - e)- 3 / 2 r 



V 



X 



o M„N U + P V P V 
2v — — V 



d) -df° + e (d T ° 



M V N V + PI (4=° - v) 



M U+1 N U+1 + P v 2 +1 



(59) 



where D.DS 6 ^ stands for the asymptotic Debye expansion up to the sixth order, we can 
write 



Ajfe, m (2,/^0;e) = Fi (z,fi; e)F 2 {z,ii,6;e) 

(consistently retaining in this product terms up to p~ 6 ). 
It is easy to see that 



(60) 



F 1 (z,fi;e)= e xp[-e-)[l + 0(e)] 



(61) 
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So, we straightforwardly get 



E £[#*,m(M,0;e = O)-F 2 (*,M;e = O)] 

n=-oo fc=0 



(62) 



oo oo 



^=00 fc=0 



^2 = Jim £ £ ex P (" e 7 ) M*, fi, 0; e = 0)(1 + 0(e)), 



(63) 



where i<2 0; e = 0). Employing the asymptotic Debye expansion for Bessel functions we 
obtain 



F 2 (z, 0; e = 0) = (2izQ n ) I + 



-3-t + 5t 2 + 3i 3 -i cos(2 0) 



+ 



+ 



6 - 6 1 + 25 1 2 + 42 1 3 + 35 1 4 - 48 1 5 - 71 1 6 
+ (-6t- 14t 2 + 12t 3 + 21 t A ) cos(20) 

-4 - 12 1 + 30 t z + 212 1 6 + 560 1 4 - 253 t° - 2044 f 



32 z/ 5 

-690 1 7 + 1562 1 8 + 852 1 9 

+ (-12 i - 84 t 2 - 76 1 3 + 406 t A + 456 1 5 - 378 1 6 - 450 t 7 ) cos(2 0) 



+ (8t 3 -3t 5 ) cos(40) 



(64) 



1 

+ - 



t 6 



- cos(0) sin(0)" 



+ 7t x 



+ 



4 ^ _ t 2)f 
t 6 r 

-8 - 20 1 + 504 1 2 + 2305 t 3 + 414 1 4 - 6696 1 5 - 17400 1 6 

128 z/ 6 L 

-14025 t 7 + 39576 i 8 + 50414 1 9 - 23856 1 10 - 32799 1 11 

+ (-8 - 168 1 - 744 t 2 + 468 1 3 + 7312 1 4 + 4842 t 5 - 16700 i 6 

-15664t 7 + 10800 1 8 + 11154 1 9 ) cos(20) 

+ (48 t 2 + 136 t 3 - 158 1 4 - 270 1 5 + 60 1 6 + 99 t 7 ) cos(4 0) 

i t 8 



— cos(0) sin(0)" 



24 + 68 1 - 91 1 2 - 169 t 3 + 63 t 4 + 100 t E 



8 Z/6(l-t2)5 

Now, equation ( |63D can be expressed in terms of the double series 



s{p, q; e) = E ^^(^ e )> *(p> ?; e ) = E ^%(^ e )> 

fc=0 fc=0 



(65) 



where 



Sp(u, e) 



-2^) E 



iQ n exp —eJ z/ 2 + f2 2 t 



v 2 + ttl 



p/2 



(66) 
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and 



T p (u,e) 



oo exp 

M £ — 



p/2 



(67) 



The series S p have been studied in Ref. ||13|| , and the series T p can be evaluated in a similar 



way. Both satisfy the same recursion relations So, they can be completely determined 
from the knowledge of iS 2 (z^, e) and % (v, e). 5 2 (za e) is evaluated in the Appendix II of Ref. 
|l3ll , and ^(v, e) can be evaluated following similar steps. 

While S p (u, e) is regular at e = for all p > 1, 7^(z/, e) is regular at e = only for p > 2. 
In both cases they are exponentially decreasing with u, making the remaining k series in 
equations ( |65D absolutely convergent, even for e = 0. 

Consequently, it is sufficient to consider the e = limit, thus getting the simplified 
expressions (see Ref. ]T3"|]) 



where 



with 



>,0) 



1 d 



[2k - 1]!! \ vdv t 



(68) 



1 



1 d 



K-l 



[2(«- 1)]!! V vdv, 



«S 2 M), 



(69) 



2// /-oo 

5i(i/,0) = — / 

7T JO 



1 _)_ e -fl/z e ^Vu 2 + l ]_ _|_ + 1 



(70) 



and 



5 2 (z/,0) 



T 2K (z/,0) 



1 _|_ e -{y,-u)/z I _|_ e (ji+v)/z 



[2(«-l)]!! (~^J T2(Z/ ' 0) ' 



(71) 



(72) 



K-l 



t ^°^[2^tf(4I;J m 



(73) 



T 2 (z/,0) 



1 



1 



1 + e~^~ v ^/ z 1 + e^ +I/ )/ 2 



(74) 



7T 



i a 



o V^ 2 + 1 \vdv 



1 



1 _)_ e -£/z e 7V« 2 +l 1 + e +v/z e ^Vu 2 +l 



, (75) 
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In particular, by inspection of the powers of p appearing in the expression of F 2 (z, p,,0;e = 
0), equation fl64j) , it can be seen that only s(p, q, e = 0) for p > 5, and t(p, q, e = 0) for p > 8 
are present in N 2 . 

These results are suitable for performing the numerical evaluations of N 2 , which will be 
described in the next section. 

Numerical evaluations become more difficult when z grows up. But, in the z > 1 region 
an asymptotic approximation valid for large z applies. 

In fact, the series in Eq. fl65|) for p > 4 can be expressed as 



(76) 
(77) 



where 



1 d 

s(p, q;e = 0) = (- 4 ^)^^^- a P-%i > 
t(p, q; e = 0) = (-Az) a M , 

oo oo 
k=0 n=0 



As discused in Ref. JIB] , a PA can be developed, for z » 1 and p > 2, as 

( 1 V 



a 



p,q 



Is 



\2nz 





r(*±i) 







1 1 • ^ 
C p — o — l, — i — 

■ H '2 2nz 



+ 



(78) 



00 (-1 



i xiH-»r(^) / i 



2nz 



r § 



2 +2 ^l 



(79) 



Notice that Eqs. ([7^,[7(|[F7]) suggest that N% is 0(,2~ 3 ) (since the terms in the series in 
Eq. (|6~2"1) are 0(1/ z) p ~ q ~ 4 with p — q = 7, the next order in the Debye expansion). So, up to 
this order, we have 



(^bag) 



87T 

IT 



t-. ~2 



/!2 + 



/i 3 \ 4// /iC(3) 



7T^ 



3 7T 90 7T 3 



(31 - 21 cos(2fl)) + 0(z~ 



(80) 



The first term in the r.h.s. of this equation is the well known contribution from the free 
fermionic field. The remaining terms come from the large z expansion of A^. For 8 = this 



reproduces the result found in Ref. [13] for the MIT-case. 



V. NUMERICAL EVALUATIONS 

In this section we will show our numerical evaluation of (N), for a selected set of values 
of z and 9, as a function of /L 

For the double series in equation (|62]) we adopt a cutoff for n and k, such that the tail 
of the series becomes negligible. Since we have taken the asymptotic expansion up to the 
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M = 6 order, the double series are strongly convergent, being the required maximum values, 
n and ko, small enough for shortening the numerical calculations. 

On the other hand, since we have analytically solved the n-sum in equation (|63|), only 
the remaining fc-sum in ( jB5|) must be numerically studied. We impose again a cutoff to these 
sums for large values of k, so as to also make negligible the tail of these series. 

In Figure [I] we show the total mean fermionic number as function of /2, for a fixed value 
of z — | and for different values of 6 < |. Similarly, Figure |2] shows the total fermionic 
number for 9 > | . 
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FIGURES 
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FIG. 1. Fermionic number for * = | — : = 0,--: = f , — : = f , --: = ^, : 
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FIG. 2. Fermionic number for z 
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The solid line in Fig. [T] corresponds to the fermionic number in the MIT bag model 
(9 = 0) |T||]. As expected, when Jjl approaches the (adimensionalized) energy of the lowest 
state of the Dirac Hamiltonian (RE = 2.04), the number jumps up to (N) = 4, which is 
the corresponding degeneracy. When 9 grows up, these four levels split into a singlet and 
a triplet state [0, the curve having a first step up to (N) = 1 at a smaller value of p,. It 
is worthwhile to remark that Nsk produces a ^-dependent (and /i-independent) shift in (N) 
(|), which is canceled by the contribution of the fermionic degrees of freedom, (Nb ag ), giving 
the expected plateau at (N) = 1. 

Finally, we show the total fermionic number for different values of z, for 9 = in Fig. |3| 
and 9 = 7r/4 in Fig. £|. These figures show how the way energy levels are populated smoothes 
when z = RT grows up . 




FIG. 3. Fermion number at 6 = 0; — :z = ^j, — -z = ^ 
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VI. CONCLUSIONS 



In previous sections we evaluated the mean fermionic number for a hybrid chiral model 
consisting in a static spherical bag containing a massless two flavored fermionic field and an 
exterior skyrmionic sector. Fermions are constrained to satisfy chiral boundary conditions 
depending on the Skyrmion profile 9(R) at the the bag surface. 

Describing the one-loop Grand partition function through the functional determinant 
of the Dirac operator subject to (chiral bag) spatial and ( ant iper iodic) temporal bound- 
ary conditions, we were able to express the mean fermionic number in the bag, (N^ ag ) = 
dG 

— -q—(/3,R,[/,;9), as a trace involving the Green function of the boundary value problem. 

This trace has been expanded as a sum over finite dimensional invariant subspaces, and 
evaluated with help of the asymptotic (Debye) expansion, allowing to show the regularity 
of the final result. 

The CCP requires to complement (N haig ) with the contribution from the tail of the 
Skyrmion, Nsk(9). This introduces a ^-dependent shift in the total mean fermionic number 
(N), which is exactly compensated by an opposite one coming from the fermionic field in 
the bag. Notice that this behavior of (A^ ag ), which reflects the spectral asymmetry induced 
by a non vanishing 9, can not be obtained by considering only a finite number of eigenvalues 
of the Dirac operator to approximate the Grand partition function. 

We would like to remark that up to now we have considered the chiral angle 9(R) as 
a free parameter. In fact, the topological charge contained in the tail of a Skyrmion-like 
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configuration, Eq. (|), depends only on 6(R) at the surface, and on the boundary condition 
linw,0(r)=O 1,0,0. 

In the description of a "baryon" ((N) = 1, since we are considering only one color) 
through this two phase hybrid model, the following picture emerges: for large bag radius 
(0(R) — > 0, i.e. Nsk — ► 0), the chemical potential \i must be fixed so that (Ab ag ) = 1. When 
R becomes smaller, the chiral angle, and consequently Nsk, rise up. The \i necessary to 
get (iVb ag )(/9, R, /i; 6) + Nsk(0) = 1 diminishes. When the radius is such that 9(R) = | the 
lowest state of the Dirac Hamiltonian dives into the Dirac sea. At this point Nsk — |, so // 
must be chosen to give (N hag ) = |. When R — ► — > 7r), the topological charge reaches 



one. Then, (A^ ag ) = 0, which is compatible with /i = 0, as can be seen from (TTSLIBLrO 



The behavior of hadronic matter at finite temperature and density has received consider- 
able attention during the last years (See for example |20[). The main motivation behind this 



effort is an attempt to understand properly the generally accepted possibility of a deconfin- 
ing phase transition from hadronic matter to the quark-gluon plasma ||21|| , with applications 
to relativistic heavy ion collisions and to the early Universe. In this framework, the study of 
the thermodynamical stability of the hybrid model considered above, through the analysis 
of the free energy, may be relevant. 

The Gibbs free energy can be obtained by integrating (iV bag ), 

G(P,R,tM;9) = F(p,R;e)- /VbagXA R, A &)dn', (81) 

J 



where F((3, R; 6) is the (fi = 0) Helmholtz free energy of the chiral bag, evaluated in ||22|| . 

When studying the system for a fixed "baryonic" number the relation between /z and 
(N) must be inverted, which requires more detailed numerical information. 

Moreover, the study of the free energy as function of j3 and R needs the knowledge of 
the detailed profile 8(f3, R). Even though one can rely on approximated explicit versions as 
in Ref. p3| , p^| , this also would require a greater numerical computational effort. We will 
report on that subject elsewhere. 
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